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ABSTRACT 



Context. Previously detected in only a few gamma-ray bursts (GRBs), X-ray flares are now observed in ~ 50% of Swift GRBs, though 
their origins remain unclear. Most flares are seen early on in the afterglow decay, while some bursts exhibit flares at late times of 10 4 
to 10 seconds, which may have implications for flare models. 

Aims. We investigate whether a sample of late time (> 1 x 10 4 s) flares are different from previous samples of early time flares, or 
whether they are merely examples on the tail of the early flare distribution. 

Methods. We examine the X-ray light curves of Swift bursts for late flares and compare the flare and underlying temporal power-law 
properties with those of early flares, and the values of these properties predicted by the blast wave model. 

Results. The burst sample shows late flare properties consistent with those of early flares, where the majority of the flares can be 
explained by either internal or external shock, though in a few cases one origin is favoured over the other. The underlying power laws 
are mainly consistent with the normal decay phases of the afterglow. 

Conclusions. If confirmed by the ever growing sample of late time flares, this would imply that, in some cases, prolonged activity out 
to a day or a restarting of the central engine is required. 
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1. Introduction 

The majority of gamma-ray bursts (GRBs) are well de- 
scribed by the blast wave model dRees & M eszaros 1992; 
Meszaros, Rees, & Wijers 19981, which details their temporal 
and spectral behaviour. In this model GRB prompt emission is 
caused by internal shocks within a collimated ultrarelativistic jet 
while afterglow emission is created by external shocks when the 
jet ploughs into the circumburst medium, causing a blast wave. 
This results in a power-law temporal decay and a non-thermal 
spectrum widely accepted to be caused by synchrotron emis- 
sion. Since the launch of the Swift satellite (Ge hrels et al. 2 004). 
it has become clear that this model for GRBs cannot, in its 
current form, explain the complexity of observed light curves 
- Swift's, fast-slew capability allows for much earlier observa- 
tions and a more elaborate picture of the evolution of the emis- 
sion, particularly in the X-ray regime using the X-ray Telescope 
(XRT, IBurrows et al. 20031) . The unexpected features detected, 
such as steep decays, plateau phases (e.g.,[Tagliaferri et al. 2005 1 
INousek etal. 20 06; O'B rien etal. 20061 ) and a large number of 
X-ray flares (e.g., Burrows et al. 2007; Chinc arini et al. 2 007 ; 
IFalcone et al. 2006] IFalcone et al. 20071 ) have revealed the com- 
plexity of these sources up to ~1 day since the initial event, 
which is yet to be fully understood. 

Prior to Swift, X-ray coverage typically began at ~ 0.5- 
1.5 days after the GRB event, and X-ray flares were detected 
in only a few bursts (e.g., GRB 970508, iPrro et al. 19981 
GRB 01 1121 & GRB 011211, IPiro et al. 20051 However, they 
are now observed in ~ 50% of all Swift GRBs, typically super- 
posed on the early light curve steep decay and plateau phases. A 



clear clustering exists in the total flare distribution in time, with 
the vast majority occurring at early times up to ~ 1,000 seconds 
dChincarini et al. 20071 IBurrows et al. 200"7l ). Various possible 
explanations have been put forward including external ori- 
gins due to patchy shells (Meszaros, Rees, & Wijers 1998; 
IKum ar & Piran 2000]), refreshed shocks (e.g., 
IRees & Meszaros 1998; Zhang & Meszaros 2002), density fluc- 
tuations (Wang & Loeb 2000; Dai & Lu 2002) or continued cen- 
tral engine activity flDai & Lu 19981 |Zhang & Meszaros 20"02] ), 
though a consensus has been slow to emerge. A very small num- 
ber of GRBs have exhibited flares on much later timescales of 
10 4 to 10 5 seconds, or approximately one day after the prompt 
event (including, in the optical, pre-Swift GRB000301C, 
Saga r et al. 2000| ). These late flares, like the early flares, are 
difficult to accommodate within the external shock model if 
the width, At, is smaller than the observing timescale, t, as is 
often the case (ILazzati et al. 2002llLazzati & Perna 20071 ). They 
can also be difficult to accommodate within the internal shock 
model because that would require prolonged activity out to late 
times or a restarting of the central engine, though a number of 
methods for doing so have been suggested (Zhang et al. 2006). 
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Here we examine the X-ray light curves of Swift GRBs for 
such late flares and compare their properties with those of early 
flares to investigate whether they are simply the tail of the early 
flare distribution or form a different sample. In section [2] we in- 
troduce our sample while in section [3] we discuss the methods 
and results of our temporal and spectral analyses. In section |4] 
we discuss these results in the overall context of the blast wave 
model of GRBs and the internal/external flare models. We sum- 
marise our findings in section|5] Throughout, we use the conven- 
tion that a power-law flux is given as F oc t~ a v~P where a is the 
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temporal decay index and f3 is the spectral index. All uncertain- 
ties are quoted at the lcr confidence level. 



2. Sample selection 

From a visual inspection of the pre-reduced Swift XRT light 
curves in the on-line repository flEvans et al. 20071 ) up to the end 
of December 2007, we identified a sample of 7 bursts which 
clearly exhibit very late (> 1 x 10 4 s) flares (FigureQ}. The def- 
inition of the cutoff at 1 x 10 4 s is entirely arbitrary but cho- 
sen so as to study a sample of the latest flares possible. As 
the redshifts of the majority of our sample are unknown, we 
are unable to correct the time of the flare to the rest frame 
time; we therefore caution that these flares are not necessar- 
ily at such late times intrinsically. However, even if we as- 
sume that these bursts are at the average Swift redshift of ~ 2.8 
dJakobsson et a l. 2006), they are still on the tail of the temporal 
distribution (Chinc arini et al. 20071) . These bursts are sampled 
well enough to allow an unambiguous identification of a flare, 
i.e., sparsely sampled bursts where the data could not rule out a 
flare were not used except for the case of GRB 0703 1 1 where the 
presence of the late flare is confirmed via optical observations. 

We only include bursts where we can unambiguously obtain 
the underlying temporal decay so that we can better constrain the 
flare parameters. This approach favours more obvious, sharper, 
stronger flares, while neglecting slower, dimmer flares. Our se- 
lection also gives a biased sample towards bursts more easily 
followed up with Swift, which are bright to late times such that 
they can be observed out to 1 x 10 5 s and a flare detected. This 
also favours bursts of shallower temporal decay as they are more 
likely to be observed out to these late times. Despite these bi- 
ases, the discovery of flares due to internal processes at late times 
would place strong constraints on models. 

The sample is summarised as follows: GRB 050502B, as 
previously discussed by IFalcone et al.l displays an energetic 
early flare as well as two overlapping late flares. GRB 050724 
displays a late flare at X-ray and optical wavelengths 
(Campana et al. 2006 IMalesani et al. 20071 respectively). We 
note also that this is considered a short burst with T90 = 
3 s over 1 5-350 keV dKrimm et al. 20051 and 2 s in softer en- 
ergy bands ( |Remillard, Swank, & Levine 2005| ). GRB 050916 
displays an obvious, sharp late time flare, though the tempo- 
ral coverage is poor. GRB 070107 displays two early flares as 
well as the late flare dStamatikos et al. 20071) . GRB 0703 11 was 
initially detected by INTEGRAL so only has XRT coverage 
past 7000 s, however, coverage at optical wavelengths has con- 
firmed the suspected flare and allowed a full study of this burst 
dGuidorzi et al. 2007atlGuidorzi et al. 2007bl ). GRB 070318 dis- 
plays both an early and late flare which was also ob- 



Table 1. The spectroscopic redshifts, z, Galactic absorption, 
A^e(Galactic) (Kal berla et al. 20051 ) and total fitted absorption, 
A^e(Total), for the bursts in our sample. 



served by UVOT in the optical (Cumming s et al. 2001) w hile 
GRB 070429A has a single late flare dCannizzo et al. 20071 ). 



3. Analysis & results 

3.1. Spectra 

The XRT event data for these bursts were initially processed 
with the FTOOL, xrtpipeline (v0 .11.4). Source and back- 
ground spectra from the Photon Counting mode data (PC; 
IHill et al. 20041 ) were extracted and bad columns corrected for, 
where necessary. Pile-up was tested for but, due to the late- 
ness of the data, was not an issue. Spectra were binned to have 
>20 photons per bin and the v(51(5 response matrices were used. 



GRB 


z 


N H (Galactic) 


N H (Total) 






xlO 22 cm" 2 


xlO 22 cm" 2 


050502B 




0.0359 


< 0.11 


050724 


0.258 1 


0.140 


75 +0 5 


050916 




0.810 


1 3 +as 
1- -0.3 


070107 




0.299 


0.47 ± 0.03 


070311 




0.236 


0.6 + 0.1 


070318 


0.84 2 


0.0144 


0.33 ± 0.03 


070429A 




0.0839 


0.22 ± 0.06 



' |Prochaskaetal. (2005| >. 2 |Chenetal. (2007) . 



The spectra, from 0.3-10.0keV, were fit with absorbed power- 
laws in Xspec 11.3.2 using x 1 statistics. The Galactic value of 
Column Density, A^h, was taken from Kalb erla et al. (2003) . 

Due to the lateness and hence dimness of these flares, spec- 
tral analysis of the individual flares similar to that done by 
Falc one et al. (2007) was not possible. Instead, we fit the spectra 
of the underlying afterglow so as to determine the electron en- 
ergy distribution indices, p, via spectral indices, /?. As the rela- 
tionships between p and ft (Zhang & Meszaios 2004) hold only 



for the underlying afterglow and not the flaring region, the spec- 
tra extracted for analysis were taken from the start of the tempo- 
ral power-law decay (i.e., after early flares or steep decay phases) 
and the late flares were eliminated. This leads, in some cases, to 
spectra with low total counts and hence poorly constrained fit 
parameters as can be seen from the results of the fits (Tables[TJ& 
0. 



3.2. Light curve modelling 

Our light curve analyses are carried out on the pre-reduced, XRT 
light curves from the on-line repository dEvans et al. 2007l) .We 
model the light curves with a combination of power-law decay 
(Table [2]) and Gaussian flares with peak, f, width, <x, and full 
width at half maximum, At = 2V21n2o- ~ 2.3548o- (Table [3). 
These fits (Figure [TJ) are used to find the relative temporal and 
flux variability (At/t and AF/F where AF is the excess flux of 
the flare over that of the underlying power-law). 

In one burst, GRB 070107, a single power-law fit was un- 
satisfactory so a smoothly broken power-law was adopted. The 
break at (1.7 ± 0.4) xlO 5 s, after which the temporal index drops 
to a = 1.9 + 0.2, does not affect our analysis as it occurs after 
the flare. In the case of GRB 050916 a FRED (fast rise exponen- 
tial decay) -like flare was a better fit than a Gaussian, returning 
a temporal power-law index consistent with that of the Gaussian 
fit. This index overestimates the final data point, a 3<x upper 
limit, but we can not imply a break from this as it is likely just 
an outlier. The FRED-like profile of the flare had very fast rise 
(~ 1800s) and decay (~ 2100s) times peaking at ~ 1.9 x 10 4 s. 
The analysis however, was carried out on the parameters of the 
Gaussian based fit to be comparable with the rest of the sample. 
GRB 070429 A, after an initial steep decay phase, has little data 
until the flare so the underlying power-law is more uncertain than 
the errors would suggest as we are unable to rule out other flares 
affecting the value. The fit also overestimates the final data point, 
which we did not include in the fit as it introduced a solution in 
the fit of the flare, not believed to give the true values of the flare 
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Fig. 1. Power-law plus Gaussian fits to the XRT light curves of each of the bursts in our sample, except for GRB 050916 which is 
better fit by a power-law plus FRED like flare. 
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parameters. The overestimate may be related to a break which 
could not be confirmed from the limited data available. 



Table 2. Indices of the underlying temporal power-law, a, the 
spectral power-law, /3, the electron energy distribution, p and the 
value of the temporal decay it predicts, a p . 



GRB 




a 


yS 






P 


a i 




050502B 


0.89 


± 0.04 


1.0 ± 


0.2 


2.0 


± 0.4 


1.0 ± 


0.2 


050724 


0.93 


±0.09 


1.4 ± 


0.5 


2.8 


± 1.0 


1.6 ± 


0.6 


050916 


0.71 


±0.05 


0.8 ± 


0.2 


1.6 


± 0.4 


0.9 ± 


0.2 


070107 


1.03 


±0.03 


1.27 ± 


0.07 


2.54 


± 0.14 


1.41 ± 


0.08 


070311 


1.2 


± 0.5 


1.2 ± 


0.2 


2.4 


± 0.4 


1.3 ± 


0.2 


070318 


1.10 


± 0.15 


1.4 ± 


0.1 


2.8 


± 0.2 


1.6 ± 


0.1 


070429A 


0.38 


± 0.05 


1.10 ± 


0.15 


2.2 


± 0.3 


1.15 ± 


0.15 



Table 3. Peak time, t, Gaussian width cr, relative temporal and 
flux variability, At/t and AF/F, of the late flares. 
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Fig. 2. Relative temporal and flux variabilities, At/t (upper) and 
AF/F (lower; Table [3]) versus the underlying temporal power- 
law decay indices (Table|2| show no correlation. 



GRB 


t 


cr At/t 


AF/F 




xl0 4 s 


xl0 4 s 


-> rv+l.i 



050724 
050916 
070107 
070311 
070318 
070429A 



7.53 d 
5.72 h 
2.040 d 
8.9 h 

15.8: 



0.18 
0.16 
0.005 
0.2 
: 1.7 



2.2 : 
1.89: 
0.12 : 

1.3 : 



0.2 

0.16 
0.01 
0.2 



0.69 : 
0.78 : 

0.14 : 
0.35 : 



0.06 
0.07 
0.01 
0.05 



17 +4 
26.27 ± 0.05 



7 +2 -° 
I 4+O.8 



1 +03 

^•^-0,9 

12 +0 - 07 



—1.9 

5 3 +L3 

-3.1 

21 +s 
70 ±40 

9+ 03 
<8.1 

1.1 ± 0.2 

1 9 +1 2 



4. Discussion 

The early behaviour of the bursts is varied; of the 5 sources with 
early observations, 2 exhibit a steep decay phase while the other 
3 display early flares over a power law decay. Since the redshifts 
of the majority of our sample are unknown, it is impossible to 
confirm the apparent agreement of the time of the flares at f avg ~ 
1 x 10 5 s, though this is certainly affected by our definition of 
late as > 1 x 10 4 s. 



4. 1 . Spectral and temporal indices and the blast wave model 

The underlying temporal indices of the bursts range from 
a ~ 0.4 to a ~ 1.2, which are quite shallow, suggest- 
ing that these flares might occur during plateau phases (e.g., 
INousek et al. 20061 IQ'Brien et al. 20061 To test this we com- 
pare the observed values of temporal indices with those derived 
from the spectral indices, assuming the standard closure rela- 
tions (Zhang & Meszaros 2004). We used these relations to cal- 
culate the electron energy distribution index, p, and the predicted 
values of the temporal slope, a p , from the measured spectral 
indices of the X-ray spectrum, (3 (Table 13. We have assumed 
that the cooling and peak frequencies are below the X-rays, 
Vc,m < vx, to estimate the shallowest slopes possible for a given 
spectral index (steeper slopes may be estimated if v m < vx < v c ). 

In all cases the predicted value of temporal decay overesti- 
mates the measured value, though in 4 of the 7 cases the values 



overlap within the lcr level, and 1 at the 2<x level. It should how- 
ever be cautioned that many of the predicted temporal decays 
have significant errors because of the low total counts of their 
X-ray spectra. In the case of GRB 070429 A, there is no overlap 
until 3.7cr, implying that it does not correspond to the regular 
decay in the blast wave model. It can be explained by energy 
injection ( No usek et al. 2 006 ) of the form E oc r 6 , which how- 
ever, would cause a break to regular behaviour which is not ob- 
served. Given the poor temporal sampling of this burst and the 
uncertainty of the underlying power-law it is difficult to make 
conclusive statements. 

In the case of GRB 070107 the overlap is similarly marginal, 
at the 3.5<x level. This is the one burst in our sample where we 
are able to say there is a break at late times (~ 1.7 x 10 5 s). This 
may be a jet break, though the observed slope of a = 1.9 + 0.2 
is shallower than the expected a p = p — 2.54 + 0.14 but may 
be rolling over to the asymptotic value expected. Another expla- 
nation may be that X-ray frequency is between that of the peak 
and cooling frequencies, v m < vx < v c , in a constant density cir- 
cumburst medium: in this case we find p = 3.54 ±0.12 implying 
a p = 1.91 + 0.06 which agrees well with our post break slope. 
The break could then be interpreted as being due to the cessation 
of continued energy injection of the form E oc f~° 6 . 

The lcr agreement of the majority of the bursts does not 
preclude the possibility of energy injection in these cases, espe- 
cially given the significant errors. It is most likely that the bursts 
are taken from both samples. The properties of the flares them- 
selves, i.e., relative flux variability, AF/F, and temporal variabil- 
ity, At/t, do not seem to be dependent on the underlying tempo- 
ral decay power-law (Figure|2]i, though with such a small sample 
this is inconclusive. 



4.2. The origin of late flares: internal vs. external shocks 



Applying 



kinematic 



arguments, 

Ioka, Kobayashi, & Zhang (2005 ) place limits on the timescale 
and flux amplitude variabilities (At/t and AF/F) allowed by 
various flare, or bump, afterglow origin models (i.e., external 
shocks): patchy shells, refreshed shocks, on-axis density fluc- 
tuations and off-axis density fluctuations of many regions. The 
limits they find are At/t > 1, At/t > 0.25, AF/F < 1. 6 At/t 
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Fig. 3. Our relative temporal and flux variabilities, At/t and 
AF/F (Table |3), of the late flares (blue error bars) plotted on the 
allowable kinematic regions (Ioka, Kobayashi, & Zhang 2005). 
Also plotted are the early (squares) and late (triangles) sample 
of |Chincarini et al. (2007| l. The lines shown are for bumps due 
to patchy shells (At/t > 1), refreshed shocks (At/t > 0.25), on- 
axis density fluctuations (AF/F < l.bAt/t) and off-axis density 
fluctuations of many regions (dashed line; AF/F < 24Af/f). 



and AF/F < 24-At/t respectively, assuming F/vF v ~ 1 and 
the fraction of cooling energy, f c ~ (v m /v c ) (/, ~ 2)/2 ~ 1/2 
( Ioka, Kobayashi, & Zhang section 3.2, section 3.3, equation 
7, equation A2). |Chincarini et al. (2007) > plot their sample's 
properties on these regions and find that while all flares may 
be explained as being due to internal shocks (At/t < 1) caused 
by long-lasting central engine activity, about half could be due 
to refreshed shocks of external shock origin and ~ 15 percent 
could only be explained by prolonged central engine activity. 

We have compared the values for At/t and AF/F of the 
late flares with those published by [Chincarini et al.l and with 
the theoretical limits of |Ioka, K obayashi, & Zhang We find that 
these properties of the late flares in our sample are in agreement 
with the distribution of values found for early and late flares 
bv lChincarini et alTI (Figure [31. In this figure, the differences be- 
tween the values for late flares in the Chincar ini et al.l sample 
(triangles) and ours (blue error bars), show that there are some 
differences between the two analyses. While the properties of 
GRB 050724 agree within our errors, those of GRB 050502B 
and GRB 050916 appear to be the most deviant. In the case 
of GRB050502B, Chincarin i et al.l fit the underlying temporal 
power-law as a broken one, while our light curve does not sup- 
port this so is fit by a single power-law. We also find that the 
early flare is better fit by two Gaussians as opposed to one. Both 
of these differences cause an offset of the properties in question. 
For GRB 050916 we find that the late flare is well fit by only one 
Gaussian as opposed to the two used bv lChincarini et alTI 

We find that all the late flares in our sample, bar that of 
GRB 070 318 with At/t > 1 an d which is also observed in the 
optical (Cummings et al. 2007) 1, may be explained by internal 
shocks, though only the FRED-like flare of GRB 050916 can be 
explained only by internal shocks. GRB 050916 lies to the top 
left of the distribution in Figure[3]and while it remains within the 
early-flare parameter space, it shows more extreme properties 
than a typical early flare. This can be seen in the shape of its flare, 
which has a very rapid rise-time and a FRED-like shape rather 



than a Gaussian. Internal shocks are almost certainly responsi- 
ble for this particular flare as density enhancements or refreshed 
shocks would not be able to create such a fast rise. The late flares 
of both GRB 50724 and GRB 0703 1 1 have also been obser ved 
in the optical dMalesani et al. 20071 IGuidorzi et al. 2007bl re- 
spectively) and this apparent achromatic nature suggests external 
origins. This is certainly consistent with their flare properties as 
shown in Figure [3] Though, these flares do lie in the expected 
range for an internal origin and are not distinct from those of 
early-time flares so this is inconclusive. 

The majority of late flares in our sample, similar to the early 
flares, may be explained both by internal shocks or one of the 
external models. The evidence from this limited sample suggests 
that these flares are no different from the sample of early flares. 
They are, most likely, late examples on the tail of the distribution 
of early flares, though this needs to be confirmed by a larger 
sample with known redshifts as it becomes available. 

5. Conclusions 

We have examined the Swift-XRT light curves of GRBs up 
to December 2007, and identified a sample of 7 bursts which 
clearly exhibit late time flares. The early behaviour of these 
bursts, where observed, is varied: either steep decay, or, flaring 
overlaid on a power-law decay. Some of the underlying power- 
law decays, at the time of the late flares, are probably due to 
continued energy injection, while the remainder are probably the 
normal decay phase of the afterglow. As the flares occur at late 
times and may obscure possible underlying breaks (energy in- 
jection, jet or spectral) we can only confirm a break in one burst. 

The evidence does not suggest that these flares are any dif- 
ferent from the sample of early flares, and hence they are most 
likely late examples on the tail of the distribution. Like the sam- 
ple of early flares, most can be explained both by internal shocks 
or external models. We should caution however, that the sam- 
ple presented here is not unbiased and there are certainly selec- 
tion effects which favour brighter bursts, more easily followed 
up with Swift, and more obvious, sharper, stronger flares. If this 
distribution is confirmed, as the sample with known redshifts and 
well sampled light curves grows, it would imply that internal 
processes produce significant flares up to a day after the prompt 
event. Hence, in some cases at least, prolonged activity out to 
late times or a restarting of the central engine is required. 
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